A diffraction-compensating 0—25 ns free space terahertz delay line for coherent 

quantum control 
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Free space delay lines provide pulses of variable time spacing for optical experiments such as 
pump-probe spectroscopy and coherent quantum control, including spin and photon echo techniques. 
However, in the terahertz (THz) region of the spectrum, beam divergence due to diffraction limits 
the useful length of traditional free space delay lines. We present a novel double-folded variable delay 
line for light in the frequency range 0.24-1.2 THz, which incorporates a symmetric arrangement of 
lenses whose spacing can be adjusted to compensate for diffraction at each delay. Scalable for use 
in other wavelength regimes, the design relays an input Gaussian beam waist to the output with 
up to 25 ns (~ 8 m) total delay and is enclosed in a desiccated volume of < 0.5 m'^. The delay 
line can deliver two or three pulses with relative amplitudes controlled via variable spacing silicon 
etalon beam splitters. Beam profiles of a 0.24 THz beam show good agreement with calculations at 
long delays, with insertion loss per delay stage of ~ 3 dB. 
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The quest for higher spectral resolution, higher sensi- 
tivity and faster time response has driven state of the 
art electron paramagnetic resonance (EPR) beyond mi- 
crowaves, to the edge of the terahertz (THz) regime^. At 
the same time, the search for a qubit with suitably low 
decoherence rates and potential scalability for quantum 
information processing has led some researchers to the 
same technologically under-serviced spectral range^''^. 
In both fields, pulsed quantum control techniques al- 
low measurements of population relaxation (via time- 
resolved decay of magnetization or emission), ensemble 
dephasing (via decay of Rabi oscillations versus pulse in- 
tensity or duration), and decoherence (via the decay of 
a spin or photon echo). Such measurements yield in- 
formation about local environments and interactions at 
the nanoscale, and are being pursued at terahertz fre- 
quencies for studies of protein conformal dynamics us- 
ing amino acid spin labels and dynamics of electrons in 
quantum confined nanostructures. We note the applica- 
tion of these techniques is predicated on the availability of 
high intensity radiation delivered in pulses of well-defined 
shape, controllable amplitude and variable time delay. 

In the optical and infrared regions of the spectrum, 
free space delay lines utilize a beamsplitter to direct a 
portion of an optical pulse along variable path length, d, 
terminated by retroreflecting mirror arrangement, which 
offsets the pulse and returns it along a parallel trajec- 
tory. The reflected portion is delayed relative to the inci- 
dent pulse by 2d/c, where c is the speed of light. At the 
other end of the spectrum, below 0.1 THz, waveguided 
microwave devices are preferred over free space delays for 
pulse generation. 

In the intermediate range 0.1-10 THz, neither ray 
optics solutions nor waveguided circuit technologies for 
pulse generation are generally applicable. This is due to 
a conspiring collaboration of failing material properties 
and the breakdown of most of the useful approximations 
that make microwave and photonics tools readily engi- 



neerable. THz sources based on semiconductor transport 
including Schottky diode multipliersi^ and quantum cas- 
cade lasersii can be modulated by microwave techniques, 
but do not yet provide adequate power for most quan- 
tum control experiments. Thus, quasi-optical sources 
and pulse generating techniques are required for most 
nonlinear and quantum control experiments^. In such 
quasioptical systems, beam diffraction and water vapor 
absorption are primary concerns. 

Here we describe a folded, dry atmosphere delay line 
which compensates for diffraction, enabling delay paths 
of up to ~ 8 m (~ 25.6 ns) for wavelengths as long as 1.3 
mm in compact 2x0.5 m footprint. The delay line is de- 
signed as dual use instrument for studies of quantum dy- 
namics of electrons confined by impurity centers^, quan- 
tum well heterostructures and quantum dots, and as part 
of a high field, 240 GHz pulsed EPR (a.k.a. ESR) spec- 
trometer under collaborative development by the UCSB 
Center for Terahertz science and Technology and the Na- 
tional High Magnetic Field Laboratory, Florida^. The 
delay line delivers two controllable amplitude pulses of 
variable delay (i.e. 7r/2, tt) for Hahn echo measurements 
(spin or photon echo) , or three pulses of equal time delay 
(7r/2, TT, 7r/2) for projecting Hahn echoes onto popula- 
tion differences, which is useful when quantum nondc- 
molition readout techniques are available for measuring 
populations directly^. 

The approach taken in the present work is reminiscent 
of an optical delay line, with the addition of a set of 
symmetrically placed lenses whose relative spacing can 
be varied to compensate for diffraction over a wide range 
of delays. The working principle is based on the fact that 
the effective focal power of a composite lens depends on 
the spacing between the elements (Fig. [1]). Consider 
a composite lens comprised of a negative lens of focal 
length — /, flanked by two positive lenses each with fo- 
cal length 2/, which are positioned equal distances away 
from the negative lens. Such an arrangement can be used 
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FIG. 1: Gaussian beam simulation illustrating the change in 
focal power with lens separation. Depicted is a system of three 
lenses of focal length 0.5 m, -0.25 m, and 0.5 m, respectively 
(curvature exaggerated for visibility). A 240 GHz Gaussian 
beam with an initial beam waist radius of 2 cm is largely 
untransformed when the spacing between the lenses is small 
(dotted lines), (b) When the spacing between the lenses is 
increased to 0.2 m (arrows), the beam (solid lines) may be 
symmetrically refocused to its original beam waist. 



to relay or "refocus" a Gaussian beam. However, when 
placed together to form a single block, the composite 
lens exhibits no significant transformation on an incident 
beam. This symmetric lens arrangement is implemented 
in our delay line because it provides diffraction compen- 
sation over the widest possible range of delays, as ex- 
plained below. 

In the present design (Fig. [2]) 45° bending mirrors 
are placed between two plano-convex 0.5 m focal length 
lenses and a —0.25 m focal length biconvex lens. The 
mirrors and negative lens are mounted to a sliding car- 
riage on a precision rail. The location of the sliding car- 
riage along the rail determines the length of the delay 
path. The sliding carriage is moved via a computer- 
controlled stepper motor and timing belt with limit and 
home switches to forestall collisions and ensure repeat- 
ably referenced positioning. The positive focal length 
lenses are mounted pair-wise on a separate computer- 
controlled slider. At each delay, the position of the pos- 
itive lens carriage can be set to relay to the output a 
Gaussian beam with a waist near the input. Other lens 
arrangements are possible for achieving the same effect, 
such as putting two or three lenses in one arm of the 
delay. However, the symmetric embodiment was chosen 
because (1) it requires only two moving stages, (2) if 
more rail sliders were needed for focusing, the additional 
thickness of the sliders would limit the minimum path 
length/delay, and (3) most importantly for long wave- 
length radiation, the symmetric design can be folded back 
and stacked on top of itself, without the addition of more 
moving stages. This is desirable to prevent a diffracting 
beam from becoming unmanageably large in the course 
of a long delay path 

The layout is simultaneously constrained by a number 
of experimental requirements. These include a desired 




FIG. 2: Schematic views of a two-level, three pulse, 
diffraction-compensating delay line. (a) 3D scaled solid 
model. The incident THz pulse (po) has a beam waist near 
the first Si etalon variable beam splitter (BSl). The first 
(undelayed) exit pulse (pi) is reflected from BSl. The trans- 
mitted portion follows a delay path of length determined by 
the location of the mirror slider (MS) along a sturdy rail. 
The lens slider (LS) is adjusted independently to refocus the 
beam such that a similar beam waist is formed near the exit. 
The return beam is split by a second beamsplitter (BS2); the 
transmitted portion exits as the second pulse (P2), while re- 
flected portion (pa) is periscoped (PRl) up to a second, iden- 
tical delay path and then periscoped (PR2) back down near 
the exit. For scale, the lenses are 100 mm in diameter, (b) 
Projected top view schematic diagram with beam paths offset 
for clarity. The dotted path represents pi, the solid path is 
P2, and dashed path is pa. The entire delay line is contained 
within a desiccated enclosure (DE). 



small physical separation of the delayed pulses (which 
limits the maximum input beam diameter) and the maxi- 
mum lens size (which simultaneously determines the min- 
imum return beam offset and minimum lens spacing) 
.These first two requirements (beam input size and max- 
imum beam diameter) fix the maximum delay that can 
be accommodated for a given wavelength. The additional 
requirements that a third pulse be twice the delay of the 
second pulse and that all outgoing pulses have similar 
beam waist diameters constrain the placement of the in- 
put beam waist location relative to the two beamsplitters. 

The optimal lens carriage position for relaying ( "focus- 
ing" ) a delayed beam is calculated using the "q" param- 
eter method^, which is valid in the paraxial approxima- 
tion. The thickness of the lenses can be treated within 
this formalism but the effect of finite lens thickness was 
found to be negligible in the present design. The cross 
section of the input beam is assumed to be a Gaussian 
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FIG. 3: Q-parameter calculation of lens separation for re- 
laying a Gaussian beam. The separation is measured as the 
distance from the center negative lens to either positive lens 
in the configuration shown in Figure [T] Lens focal lengths 
are 0.5 m and -0.25 m. Minimum beam waists are chosen for 
each frequency to yield real solutions over the entire range of 
delay. 



intensity profile, and is characterized by two parameters: 
beam waist radius, and radius of curvature. The "focal" 
condition requires the outgoing, delayed beam have the 
same beam waist radius and radius of curvature as the in- 
cident beam. This is sufRcient to determine the required 
spacing between the positive and negative lenses. The 
focal lengths of the lenses are chosen such that the spac- 
ing required to focus at short delays is not less than the 
minimum allowed by the interpositioned bending mirrors. 
The required lens separation is then calculated given the 
delay, the lens focal lengths, and minimum beam waist 
radius (Fig. [3]). 

The maximum delay per delay line is 12.8 ns (25.6 ns 
total). The minimum delay per delay line is limited by 
the minimum lens spacing, and is dependent on operat- 
ing wavelength, falling in the range 3.5-4 ns. The lenses 
are kinematically located and can be removed from the 
beam path for visible alignment, and for use at frequen- 
cies which do not require diffraction compensation. The 
delay line contains two variable beam splitters (Fig. [51 
BSl and BS2), which can be removed or replaced with 
fixed ratio wire mesh beam splitters or fully reflecting 
mirrors. A variety of configurations are possible. With 
the second beamsplitter removed, the delay line produces 
two pulses of relative delay 3.5-12.8 ns. With a mirror 
at the position of the beamsplitter, two pulses are de- 
livered with delays of 7-25.6 ns. Additionally, an extra 
3.5 ns compensating path was designed with diffraction 
compensating lenses, which can be inserted in the path 
of the undelayed beam (Fig [21 Pi ) to produce pulse of 
delay 0-9.3 ns or 3.5-22.1 ns. With the variable second 
beamsplitter in place, three pulses are generated, which 
are separated by equal, variable delays of 3.5-12.8 ns. 

For pulses of nanosecond duration, an etalon-based 
variable beam splitter is used to control the relative inten- 
sities reflected by the first beamsplitter (Fig. [21 pi) and 
transmitted by the second beamsplitter (Fig. [H ps). The 
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FIG. 4: Si wafer etalon transmission at 240 GHz. The solid 
line represents the theoretical transmission of two 500 fim 
thick Si wafers at 45° angle of incidence, the open circles are 
experimental measurements of the beam power normalized to 
the power with the wafers removed. The only fit parameter 
is the horizontal axis offset. 



etalon is comprised of two closely spaced 500 ^m thick 
Si wafers at 45 degree incidence. The wafer spacing can 
be adjusted via a translation stage to achieve transmis- 
sion in the range 0.1-1 (reflectivity 0-0.9) (Fig. [J). For 
shorter pulses in the tens to hundreds of picoseconds time 
scale, fixed beam splitters based on wire meshes should 
be used, as multiple refiections from the four interfaces 
in the etalon would distort and stretch the pulse. 

The lenses are cut from 100 mm diameter polypropy- 
lene (PPE) round stock on a CNC lathe, and the mirrors 
are fabricated in a cost-effective manner by e-beam de- 
positing 600 nm of Au on 650 ^m thick polished Si wafers 
with a 20 nm Ti or Cr adhesion layer. During operation 
at frequencies above 0.5 THz, a diaphragm pump is used 
to continuously circulate the enclosed 0.5 m^ volume of 
air through an air dryer/filter, such as is commonly used 
in compressed air lines. Atmospheric moisture, as mea- 
sured by a NIST traceable hygrometer, is reduced from 
ambient 40% to less than 1% relative humidity in approx- 
imately 45 minutes. 

The delay line is tested with 240 GHz solid state fre- 
quency source based on Schottky diode multipliers that is 
coupled to a free space TEMOO Gaussian mode via a cor- 
rugated feedhorn. The source passes through two lenses 
which match the Gaussian beam to produce a minimum 
beam waist diameter of 2 cm at the input of the de- 
lay line (Fig. [2l BSl), which enables us to examine the 
effect of unmodeled nonidealities. The most significant 
nonidealities include (1) alignment error, (2) input beam 
parameter uncertainty, (3) wavefront distortion due to 
thickness dependent loss in the PPE lenses, (4) succes- 
sive clipping of the beam peripheral fields by mirrors and 
lenses of size on the order of the beam waist diameter, 
(5) the paraxial approximation, and (6) lens and mir- 
ror defects due to material inhomogeneity, manufactur- 
ing defects and strained mounting. The delayed beam 
is profiled via a rastering pyroelectric detector of size 2 
mm on a side. The metal housing is carefully shielded 
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FIG. 5: Delay line beam profiles at 240 GHz. Beam pro- 
files are obtained at (a) input, (b) p2 and (c) pa outputs 
as a function of delay. Profiles are measured by rastering 
a pyroelectric detector of size 2x2 mm over a 5.2x5.2 cm 
area. The slight vertical elongation apparent at long de- 
lay settings (e.g. (c)) is due to a known mirror imperfec- 
tion. Profiles are shown in normalized linear grayscale for 
visibility, (d) For comparison, centered horizontal and ver- 
tical cross sections of figure [5j; are shown with correspond- 
ing traces from a 2D Gaussian least squares fit of the form 
a\ exp — [((a; — a;o)/(ia;)^ -f ((j/ — yo)/rfy)^], where dx= 2.3 cm 
and dy= 2.4 cm. 



with black THz absorbing foam to prevent feedback on 
the source. Prior to passing through two input match- 
ing lenses, the beam is indistinguishable from a perfect 
Gaussian. Afterward, the center of the beam profile is 
slightly flattened due to loss from passing through sev- 
eral millimeter thick 0.25 m and -0.5 m focal length PPE 
lenses (Fig. [5^). 

Figure[5)3-c show the output beam profiles at locations 
P2 and p3 in Figure [21 at short and long delays, respec- 
tively. The beam width of a series of beam profiles, es- 
timated by a fit to a Gaussian profile (e.g. Fig. 
and plotted in Fig. [6^, alongside the integrate intensity 
(Fig. [Hb). At short delays over focusing is noticeable, 
relative to the expected result from q-parameter simula- 
tions. This is attributed to distortion of the beam inten- 
sity profile by thickness-dependent loss in the polypropy- 
lene lenses. The difference in lens thickness between the 
center and edge of the beam is most pronounced in the 
positive lenses, where the beam is large. The center part 
of the beam experiences more absorption, which results in 
a flattened intensity proflle that requires less diffraction 
compensation than the Gaussian proflle assumed by the 
Q-parameter method. The effect is only partially com- 
pensated by the negative lens, since the beam diameter is 
much smaller at the negative lens. At longer delays, the 
return beam is slightly wider than expected. This is due 
to the far held diffraction from non-Gaussian components 
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FIG. 6: Delayed beam properties, (a) Average exit beam 
waist diameters (dx + dy, from 2D Gaussian fit explained in 
figure [5)1 caption) of a series of beam profiles are plotted ver- 
sus delay. Measurement error from repeated measurements fit 
within the respective labels. Asterisks label beam waists mea- 
sured with the lens separation decreased relative to the ideal 
calculated position of figure O (b) Transmitted power, mea- 
sured as the integrated beam profile intensity, is normalized 
to the power in the p2 output of the delay line at minimum 
delay and plotted versus delay. The difference between ps 
and p2 at equivalent delays refiects the insertion loss (~ 3 
dB) of the second delay line (including three lenses and seven 
mirrors). 



(i.e. > 1) in the intensity proflle introduced by the 
thickness dependent lens loss, as well as clipping of pe- 
ripheral flelds by mirrors of size on the order of the beam 
diameter, as evidenced by the negative slope in figure [6^ 
at delay times corresponding to long delay paths. 

The total power in the delayed beams is reduced by 
a factor of '^0.5 per delay stage, attributable primarily 
to reflection from the six lens surfaces which are not an- 
tireflection coated, secondarily to lens absorption^*', and 
lastly to diffraction losses. While this seems large in com- 
parison to losses achievable in the optical regime, the util- 
ity of a compact, desiccated delay line becomes apparent 
when compared to the typically much higher losses of 
THz beams traversing the similar distances in a nondcs- 
iccated atmosphere. 

The delay of short pulses is demonstrated via use of 
fast PIN diode to switch the 15 GHz microwave oscilla- 
tor that seeds the multiplier chain of the 240 GHz source, 
producing single triggered pulses of 10 ns duration and 
~10 mW average power, which are detected via a Schot- 
tky diode (Fig. [T]). 

Future quantum control experiments will employ the 
existing UGSB free electron laser (FEL) as a frequency- 
tunable source providing kilowatts of peak power, to- 
gether with an optical pulse slicing system — which deliv- 
ers single THz pulses of variable (3 ps-3 ns) duration at a 
1-10 Hz repetition rate. The sliced FEL pulse will enter 
the delay line and emerge as two or three THz pulses of 
similar size, with variable intensity and delay suitable for 
Hahn echo and other pulsed coherent control techniques. 

We emphasize that the delay line design presented 
herein is optimal for a subset of applications requiring 
time delayed THz pulses. There exist a variety of tech- 
niques for THz pulse generation and delay based on the 
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FIG. 7: Temporal output of delay line at 25.6 ns delay. A 10 
ns duration 240 GHz pulse enters the delay line at po in figure 
[2] and results in three output pulses (Fig. [H pi, P2 and ps). 
The three output pulses are directed onto a single Schottky 
diode detector and the temporal response is recorded by an 
oscilloscope with 1 GHz limiting bandwidth. The relative am- 
plitudes of the pulses are controlled by variable beam splitters 
(see Fig. illustrating the wide dynamic range of reflection 
and transmission that can be achieved. 

timing of optical pulses, including photo-activated semi- 
conductor switchesiS'i'^, and THz generation via ultrafast 
pulses in nonlinear crystals^^. When possible, THz pulse 
generation based on optical pulse timing is preferred be- 
cause optical delay lines do not required diffraction com- 
pensation. In situations where optically gated technolo- 
gies cannot produce the correct pulse shape or timing, a 
number of methods for physically delaying THz beams 
can be used. For fixed delay applications, diffraction in 
a THz beam can be managed via waveguides based on 
polymer, metal or semiconductor wires or ribbon a^^i^^'^^ . 
For power hungry applications, intolerable waveguide di- 



electric losses and insertion loss are the primary concern. 
In long travel Michelson interferometers, arrays of corner 
cube retroreflectors compensate for diffraction by piece- 
wise retroreflecting the face of a divergent beam. How- 
ever, Michelson interferometers back reflect half the in- 
cident power and provide no offset of the return beam, 
resulting in a train of decaying pulses, rather than di- 
viding the power between only two pulses. Deforniable 
mirrors are another technology which may potentially be 
leveraged to provide diffraction compensation. 

Future implementations of a diffraction compensating 
THz delay line will benefit by using Fresnel lenses, which 
can be thinner, and will cause less distortion of the Gaus- 
sian beam profile from thickness-dependent loss. Addi- 
tionally, for nonideal beams where diffraction cannot be 
satisfactorily managed to produce a constant intensity of 
pulses versus delay, a method for compensating a change 
in power with changes in delay should be implemented. 
One such option would be inserting a variable attenu- 
ator/reflector in each of the two delay lines. Another 
alternative is to use one variable attenuator and motor- 
ize the etalon beamsplitters to compensate for measured 
changes in output power. An alternative layout suitable 
for 2D NMR-type experiments in which the delay of a 
third pulse is varied independently of the second pulse, 
could oppose or juxtapose two independent delay lines, 
with similar beam splitters and sliding lens carriages. 
These solutions coupled with the presented method for 
diffraction compensation will provide suitable pulses and 
delays for coherent control experiments in the sub-mm 
and THz spectral ranges. 
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